Canonical Strangeness and Distillation Effects 
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Strangeness canonical ensemble for Maxwell-Boltzmann statistics is reconsidered 
for excited nuclear systems with non-vanishing net strangeness. A new recurrence 
relation method is applied to find the partition function. The method is first gen- 
eralized to the case of quantum strangeness canonical ensemble. Uncertainties in 



^ . calculation of the K + /ir + excitation function are discussed. A new scenario based 

I/-) | on the strangeness distillation effect is put forward for a possible explanation of 



anomalous strangeness production observed at the bombarding energy ~ 30 GeV. 
The peaked maximum in the K + /tt + ratio is considered as a sign of the critical end- 
point reached in evolution of the system rather than a latent heat jump emerging 
from the onset of the first order deconfinement phase transition. 

I. INTRODUCTION 



$_i ■ In recent years, statistical thermal models have widely been used for analyzing heavy 

<s ■ n 

ion collisions (see the review article These models describe available phase space at 
the final stage of nuclear collision when inelastic interactions between hadrons cease, i.e., at 
the chemical freeze- out. The grand canonical (GC) ensemble is defined in the large volume 
limit, so the volume parameter V is an extensive coefficient in the expressions for such 
quantities as mean particle numbers controlled by the chemical potentials. Application of 
the GC ensemble also implies that the number of produced particles carrying a conserved 
charge is sufficiently large. In this description the net value of a given charge (e.g. electric 
charge, baryon number, strangeness, charm etc.) fluctuates from event to event. These 
multiplicity fluctuations with respect to mean particle multiplicity can be neglected only 
if the particles carrying the charges are abundant, i.e., for the system of a large volume 
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and/or high temperature. Here the charge is conserved on the average and, therefore, the 
GC description is adequate. 

In the opposite limit of low production rate, fluctuations of the particle number can be as 
large as its event averaged value. In this case, the charge conservation has to be implemented 
exactly in each event. So the ensemble to be chosen is the canonical one. Then, the V- 
independent fugacities are replaced by nonlinearly V-dependent canonical chemical factors. 
The exact conservation of quantum numbers introduces certain constraints on the system 
considered. In the formulation of thermodynamics of strongly interacting matter one needs 
to implement these constraints that are originally governed by internal symmetry of the 
Hamiltonian 0,0,0]. Consequently, the equilibrium distribution of particle multiplicities 
can differ from that expected in the GC limit. 

In principle, for an equilibrium system all conserved charges should be treated canoni- 
cally. However, in a good approximation, baryonic and electric charge conservation can be 
considered grand-canonically even for high-energy pv collisions. This treatment results only 
in a slight overestimation of the canonical effects p|. On the other hand, the importance 
of exact strangeness conservation for hadron production in heavy ion collisions has been 
clearly demonstrated some time ago in Refs. jfj Q], where closed expressions for suppression 
of chemical factors were derived in terms of classical Maxwell-Boltzmann statistics for the 
global strangeness neutrality condition S — 0. 

In this paper, we reconsider the derivation of the strangeness constraint for nuclear sys- 
tems with non-vanishing net strangeness S applying two different methods for evaluating the 



partition 
relations 



imction and ensemble av erag es: the direct method and the method of recurrence 
In the line of Refs. P, licl |. the second method allows us first to generalize our 
canonical strangeness evaluation to the case of quantum statistics. Two Appendixes con- 
tain details on the calculation of canonical quantities by recurrence relations. The derived 
formulae will be used to analyze the kaon excitation function from heavy ion collisions and 
particular attention will be paid to the discussion of an anomalous maximum of excitation 
functions at the beam energy near 30 GeV. A new possible scenario of its formation as a 
manifestation of the critical end-point is advanced. 
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II. CLASSICAL STATISTICS 

Let us construct the thermodynamic partition function Zs with the exact strangeness 
conservation for a relativistic perfect gas of N hadron species at temperature T in space 
volume V. Every hadron species a (a = 1,...,N) is characterized by the baryonic b a , 
electric q a and strangeness s a charges, the mass m a and the spin degeneration factor g a . 
Let us assume that this particle ensemble is described by the Maxwell-Boltzmann statistics. 
In the occupation number {is ap } representation, the constraint for an exact conservation 
of the net strangeness S can be directly imposed on the grand-canonical partition function 
conserving the average baryonic and electric charges : 

Zs = £ 6(£ s a v ap - S) — ^ e-^w* , (1) 

a* 11 "ap- 

with 

£ ap = \fp* + m 2 a - [i B b a - [i Q q a . (2) 

Here (3 = 1/T, //B,/ig are baryonic and electric chemical potentials, respectively, and multi- 
index p = (p,a,.. .) includes particle momentum p, spin projection a and other quantum 
numbers. The summation in $Q) is carrying out over all permutations {v a p} of occupation 
numbers. In this representation, the ensemble average for the operator 0({v ap }) is defined 
as follows : 

(0}s = ±- £ *(E s ^p - S ) Tf^ e^ E -^ £ap ■ (3) 

These basic expressions will be used below for finding the mean occupation numbers and 
mean particle multiplicities. 

A. The direct method for solving the partition function 

The partition function Zs can be summed up over the momentum p by inserting the 
identical unit in the form J2{n a } Ua 3{J2 P v ap — n a) = 1 into Eq.fJH). Then, the partition 
function is reduced to 

z s = £*(£*«*« -s) n^r ' ^ 

{n a } a a n *- 
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where the grand-canonical single particle partition function 



-Pe 



l3{pB b a +fiQ g a ) 



(2tt) 3 



d 3 p e ~^^ +ml 



(5) 



Using the integral representation for the Kronecker symbol in Eq-Q, one can carry out 
summation over all particle numbers {n a } by means of series expansion of the exponent. So 
we have 

1 2 



2tt 



cxp 



iScj) + £ Zle 



is a (p 



(6) 



One should note that this representation of Zs based on the classical statistical partition 
function (|T} coincides exactly with the starting expression used in Q, 0, Q]. 

Let us proceed now from the summation over species index a to the sum over all available 
strangeness values s for individual particles in the system. So we substitute into Eq.fJHJ) the 
identical unit in the form J2t= a - Smax ^s,s a = 1, which is valid for any species a under condition 

I &oi I — Smax 



b 2tt 



2tt 



exp 



is<fi+ Y s * 



(7) 



where 



S s — b s ^ Sa Z) y 



S r 



< s < s r 



(8) 



By rewriting the exponent of the sum over — s max < s < s max in Eq. f T j l as the product of 
exponents and expanding every exponent into series, e x = ^x n /n\ t we get 



exp 



Y ex p 

{Ns} 



Y^n 



n 



NA 



(9) 



where N s are non-negative integer numbers. Here the following short-hand notation has 
been used for the sum Z){7v s } = X)jv_ s ■ ■ -J2N Smax - After substitution of Eq.^ into JH) 
and integration with respect to 0, the strangeness canonical partition function Zs is reduced 
to the following equation : 

Smax Smax §~^ s 



Zs= £<*( £ sN.-S) n 

■[7V S } ® — Smax S — S m 



NA 



(10) 



If the integral representation for the Kronecker symbol is used in Eq. fflO^ . we come back 
to Eq.0 for the partition function in the statistical model 
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A number of sums over the set {N s } in Eq. (fT0|) can be reduced. For this one should 
introduce the identity 

Smax OO 

17 E S(N s -N„ s -n s ) = l (11) 

s=l n s =—oo 

into (fTU|) . to be valid for every fixed number set {N s }, and change the order of summing with 
respect to {n s } and {N s }. Then, using the Kronecker symbol in Eq. (fTT|l . the summation 
for the set {N s } can be done explicitly. It is noticed that these sums are grouped by pairs 
for every =f|s| and every pair results in the modified Bessel functions. Finally, for Eq. lfTUl) 

we get 

CO OO ^max $max 

Z s = e So E ... E 6(^sn s -S) n< s W* s ), (12) 

tl\ — — oo n Smax =—co s=l s=l 

where /„ is the modified Bessel function and variables 



■1 s 



S-s S s , a s — \J S s / S_ s , s — 1, . . . , s max . (13) 



Taking into account Eq.((HI) we see that x s ~ V while a s are independent of the system 
volume. In specific calculations the maximal hadron strangeness s max is not large, which 
simplifies essentially Eq. (fT2|) . since the number of infinite sums is (s max — 1). In particular, 
for s max = 1,2,3 the partition function (fT2|) can be expressed as follows : 

Z s = e 



So E E «f" 2n_3m a n 2 a™ I s -2n-3m(xi) I n (x 2 ) I m (x 3 ) , s max = 3 ,(U) 



n=—oo m=— oo 
oo 



Z s = e So E a l~ 2n «2 h-2n(xi) In(x 2 ) , s max = 2 , (15) 

n=— oo 

Z s = e So af IsM , s max = 1 , (16) 

where the symmetry property of the modified Bessel function i_ n = I n has been used. 
Associative production of multi-strange baryons is neglected in Eq. lfTBl . This simplified 
equation for S ^ in the case of s max = 1 has been derived earlier in |l|. The result for 
W = 4 can be found in Q. The case of fol. caaonica, treatn.at of beta baeyoa and 
strangeness S ^ in a little bit different technique was given in [12]. However, the electric 
charge conservation is taken into account by neither canonically, nor grand canonically. 
It simplifies calculations but the final result differs from our Eqs. (fT^|) - (fT6|) . The single 
particle partition function Z\ given by Eq.(P can be integrated resulting in the well-known 
expression 

Z\ = K 2 (Pm a ) eflMB^+waO (i 7 ) 

2tt z p 

with the Bessel function K 2 . 



6 



B. Method of recurrence relations 

Let us subdivide the set M of particles of all species into three sub-sets M± and M 
according to the value of individual strangeness charge s a of involved hadrons. For a fixed 
set of occupation numbers {v ap } we introduce notation of positive S + and negative S_ 
components of the net strangeness S of the system 

E E 1 s «\ u <*p = s t , ( 18 ) 

aeM T P 

where M + and M_ are the sets of all species a for which s a > and s a < 0, respectively. 
For the zero component the following equality is fulfilled Y,aeM J2 P \ s a \vap = 0, since M is 
the set of all non-strange hadrons, s a = 0. Let us substitute into Eq. (fTJ) the product of two 
identical unit operators corresponding to a fixed set of occupation numbers {z/ Qp } : 

oo 

E 5 (E £M*v-&f) = i. (19) 

5 T =o oeM T p 

Then for the strangeness canonical partition function we have 

oo oo 

Zs = £ E 8(-S- + S + -S) z£> 4: } , (20) 

S+=0 5_=0 

where the partition functions for positive Zg~) and negative Zg_ components of the net 
strangeness S have the following form in the case of classical statistics : 

4? = E S( £ £ \*a\ »« P - S T ) = \- , e~^^ X> "*> . (21) 

{u ap UeM T P J^p aP ' 

The positive S+ and negative S- strangeness components take integer non-negative values. 
The appropriating partition functions can be represented as follows (see Appendix A) : 

Z S T = "FT" E \ S a\ Z a Z s \ Sa \ = E sS Ts Z S T -s , (22) 

where the second equality in (j2*2*j) was obtained by introducing once more identical unit in 
the form J2t=T $s,\s a \ = 1, to be valid for any species a under condition \s a \ < s max . Here 
by definition = 1, s max = max{|s a |} is maximal modulus value of particle strangeness 
in the given system and auxiliary functions are defined as 

S TS = E ^,\s a \ Z x a , 1 < s < s max . (23) 
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It is of great interest that mathematical structure of Eqs.(J22|) allows one an exact recursive 
solution because the partition function of the given strangeness component Zg is expressed 
through the sum of partition functions but with smaller index of the strangeness compo- 
nents. This equation can be solved iteratively starting from the lowest component Zq = 1. 
For calculating Z<^ this procedure converges quickly as is seen from Eq. ([2*T| : High values 
of S T in the argument of the delta function should be compensated by large occupation 
numbers which are exponentially suppressed with increasing negative exponent. This sit- 



uation reminds very much of 
multifragmentation model 



of find ing the exactly solvable partition function in the i 
P 111 13 and in the context of quark-gluon plasma 



"ind ing the exactly solvable partition function in the nuclear 

where 

the method of recurrence relations turned out to be very powerful. 

The partition function Z^°> with zero component of net strangeness S is reduced to the 
following generating function Z^> = e s ° with Sq = Y^aeMo according to Eq. (|23|) . 

C. Ensemble averages 

The mean occupation numbers can be found as ensemble averages (jBJ) for the operator 
0({v a p}) — v ap . Using mathematical technique for summing over occupation numbers 
as developed in 0, 1^|, we get for the mean occupation numbers the following result (see 
Appendix A) : 

(*V>s = %^ . (24) 

The mean particle number of species a is obtained by summing the mean occupation numbers 
(JH| over momentum p 

Ms = ZMs = Zl ^ • (25) 

Eq. (|25|l was given also in [11]. It is easy to convince ourselves that the mean number of 
particles with strangeness s in the system with the net strangeness S can be represented as 
follow : 

(N s ) s = £ 8 s>Sa (n a ) s = S S ^ (26) 

with S s defined by Eq. (|2Hj) . As a particular case, for the particle density in strangeness 
canonical ensemble n a = (n a )s/V we have: 
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n a = -f— X! E ^ s ^ 2n " 3m a2«™^s Q -2n-3 m (xi)/ n (x 2 )/ m (x 3 ),W = 3,(27) 

S n=— 00 m=— 00 
7l e 5o 00 

n„ = -f— at Sa ' 2n a n 2 I s ^ 2n ( Xl ) I n {x 2 ) , s max = 2 , (28) 
V Zjq 

° n=—oo 



_ z\ ( y j s -. g (gi) _ , . 

Similarly to Eqs. fTlll - fTHll only Eq. (|2Tjl takes into account that strangeness of, for example, 
kaons with s a = 1 can be compensated associatively by production of any multi-strange 
hyperon. For non-strange hadrons, s a = 0, the particle density is 

n a = Q (30) 

as follows from Eq. ([25|) . For vanishing net strangeness, 5 = 0, Eqs. ([27]) - ([29|) coincide with 
those in Ref. |l] (see also references therein). 



III. QUANTUM STATISTICS 

As compared to Eq.JTJ, in quantum statistics the partition function with exact conser- 
vation of net strangeness and with in average conservation of baryonic and electric charges 
can represented as follows : 



Zs = E 5 (E s ^a P ~ S) e"^ 1 ^ £ap , (31) 

where e ap is defined by Eq.(j2J) and occupation numbers take values u ap = 0, 1, . . . , K a with 
K a = 1 for Fermi-Dirac statistics and K a = 00 for Bose-Einstein statistics. By means of 
Eq.fjnj]) the partition function (|3Tf is reduced to Eq. (|20| with partition functions for positive 
and negative strangeness components given by equation 

4? = E 8( E E KK P - S T ) e~^^ ^ e - . (32) 

As is shown in Appendix B, the following recurrence relations are fulfilled for these partition 
functions : 

^ ^ J [%/l*a|] I ) 1 SmaX ' 5T//S ' ( ) 

= 7^ E E l Sa I Z a Z S T -l\s a \ = 7T~ E E S $Ts,l Z S^-ls > (33) 
°T tt gM T 1=1 °=F s=l i=l 
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where = 1 and new auxiliary functions 



Sts,i = E 5 s, M Zl (34) 



=F 



Z l a = ^E e_/3kQP • ( 35 ) 
p 

The quantity y a i = (— for Fermi-Dirac statistics and y a i = 1 for Bose-Einstein statistics. 
Note that in the case of Maxwell-Boltzmann statistics y a i = 5^\. Eq. (|3~3l exactly follows 
the structure of the corresponding Eq. ip^l for classical statistics and can be solved by the 
recurrence relation method. 

The partition function for zero strangeness component looks like (see Appendix B) 



Z^=exp(Y,So,i) , (36) 
i=i 

with 5*0,; = J2a&M Z l a . Notice that under condition |e _/feap | < 1 the partition function 
can be represented in more familiar form : 

Z<® = II IK 1 ± e-^f 1 , (37) 

a£M P 

where the upper sign corresponds to Fermi-Dirac statistics and the lower sign does to Bose- 
Einstein one. 

As shown in Appendix B, the mean occupation numbers for hadrons with negative and 
positive strangeness, a G M T , are given as 

7(0) oo oo [%/|s Q |] 

Ms = ~y- E E K-S- + S + - S) 4? E y a ie-^zgl lM ,aeM T ,(38) 
^s s+=os =0 i =1 

while those for non-strange particles are 

oo 

(v ap )s = e ~ Pl£ap > « G M . (39) 

i=i 

Under the same condition |e _/3eap | < 1 we arrive at the familiar expression for mean occu- 
pation numbers (|3~9*j) : 

' V «P/S = e (3e ap ± i 

To get the average number of particles a one should sum the mean occupation numbers 
(|38|) and (|39|) over momentum p. For strange particles a G M T we have 

7(0) oo oo [S T /|s Q |] 

Ms = ~y- E E 8{-S- + S+-S) Zf2 E Z i 4?-i[.al > a G M T . (41) 
^s s+=os _ =0 i=1 



Ms = Hf ^ 1 , a G M . (40) 
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For non-strange hadrons a G M the following relation is valid : 

oo 

(n Q } s = XX, aeM . (42) 
1=1 

After integration the single particle partition function (|35|) is reduced to (cf. with Eq. (fT7|) ) 

Z l a = y a i^jf K 2 (pim a ) e ^» b "+^ . (43) 

One should note that all formulae for Maxwell-Boltzmann statistics from the proceeding 
section are recovered by the simple substitution y a \ = S^i into appropriate quantum rela- 
tions. 

IV. STRANGENESS EXCITATION FUNCTION 

To reveal canonical effects in the obtained statistical relations, let us consider excitation 
functions for strange particle multiplicities. These multiplicities are calculated along the 
"experimental" freeze-out curve defined at any bombarding energy by the temperature T 
and the baryon chemical potential /jb '■ 

1270 170 

H B [in MeV] = 1 + T[in MeV] = - - e _ , 48( ^_ 3 8) , [v^ is in GeV], 



as approximated in [l|,lly|. The perfect hadronic gas includes all mesons with m a < 1.6 GeV 
and baryons with m a < 2 GeV. In calculations of a particular hadron yield the contribution 
from decays of heavier hadrons is taken into account in addition to the primary yield 

(n a ) = (n a y rmar y + ]T fir ( 7 - a)(n y ), 

7 

where the branching ratios Br{p/ — > a) are taken from the Review of Particle Physics 
Certainly the final answer will depend on whether only strong decays or also weak ones are 
considered. The latter circumstance is controlled by the experimental acceptance 

The K + /tt + excitation function is shown in Fig.l for systems formed in high-energy 
heavy-ion collisions. The curve for the strangeness neutral system described by the classi- 
cal Maxwell-Boltzmann statistics and without taking into account weak decays reproduces 
closely the earlier statistical model results well as the review 111). Note that 

two curves calculated within the classical ensemble by two above-discussed methods coin- 
cide exactly and are indistinguishable in the figure. Calculations have been carried out using 



11 



a complete equation pTj). If s max = 1 is assumed (i.e., Eq. (|29|) is applied), which physi- 
cally implies that associative production of the K + mesons is accompanied by creation of 
single-strange-charge particles only, the K + /n + ratio decreases by about 5 % at the colliding 
energy y/s ~ 8 GeV. 



0,30 



0,25 - 




s 1/2 [GeV] 



FIG. 1: Excitation function for relative strange particle production, K + /ir + , in heavy ion collisions. 
The symbols are experimental An data of E866 (triangles), NA49 (squares) and BRAHMS (circles) 
collaborations Q> Q|- The lines are thermal model calculations for two types of statistics: 
classical (solid) and quantum (dotted) for the net strangeness of the system S = and classical 
Maxwell-Boltzmann ensemble for systems with non- vanishing net strangeness S = 5, 10 (long 
dashed, dash dotted) and S = -10 (short dash dotted), respectively. The radius of the strangeness 
correlated sphere is R = 7 fm and electrical chemical potential /j>q = and s max = 3. The additional 
contribution into K + and ir + from strong decays of hadrons has been taken into account. 



The inclusion of quantum statistics slightly decreases the K + /n + ratio at > 6 GeV. 
As is found, the contribution of weak decays of hadrons with the life-time r < 1CT 10 sec 
decreases this ratio by (15)% without changing the shape of the excitation function in 
agreement with recent finding in 
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It is noteworthy that due to local strangeness conservation, two volumes are usually 
introduced. The first V is coming as a normalization factor of the partition function (see, 
for example, Eq. (f30| . which is dropped out if the particle ratio is considered. The second 
volume, V c , enters the argument x of the modified Bessel functions I n (x)in Eqs. (|2Tj) - (|2H 
and characterizes the range of strangeness correlations. It was shown |l| that for heavy 
ion collisions V c is proportional to a number of participant nucleons. The presented results 
correspond to a constant strangeness correlation volume V c treated as a free parameters. 
Following Ref.j^l, this volume is taken as a sphere with the radius of about R ~ 7 fm, 
which ensures the agreement with experiment for the fast growth of the K + /tt + ratio at 
moderate energies. Asymptotically, this ratio goes to the GC results being a constant (or 
slightly depending on a particular choice of the freeze-out state). However, in dynamical 
models this fall may be energy dependent. In particular, in the expanding fireball model (21| 
the strangeness correlation volume was associated with that of the initial thermalized state 
as the Lorentz-contracted cylinder V c = ttR 2 ■ 2R/ r y cm , where ■y cm is the 7-factor of colliding 
nuclei in the cm. system. In this case, for high colliding energies the K + /it + excitation 
function falls down faster than for a constant volume. Nevertheless, the observed peaked 
maximum in the excitation function is not reproduced within this dynamical model (see 
also the comprehensive analysis It is remarkable that such a behavior has really 

been predicted in terms of a simple statistical model of the early stage as the onset 
of the first order deconfinement phase transition. It is of interest to note that the above 
mentioned calculations, made in a similar scenario but in a more elaborated model of an 
expanding fireball 2l| with various equations of state, do not give rise to a peaked maximum 
and calculated excitation functions are practically independent of the equation of state. In 
other words, the condition of global strangeness neutrality at the freeze-out washes out 
particularities of strangeness evolution at the earlier stage. 

Generally speaking, the introduction of two volumes, V and V c , into the statistical de- 
scription based on the empirical freeze-out T—\xb curve breaks thermodynamical consistency 
of the model. Indeed, in the presented calculations there is only V. On the other hand, the 
T and /zg values were extracted from the particle ratios and the volume is not fixed. To 
correct the description we introduce explicitly the conservation of the total baryon B and 
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electric Q charges 

B = EK> & - ( 44 ) 

a 

Q = J2Mla ■ (45) 
a 

Since //# is known from the T — [Ib curve, the first equation ijHJ) defines the volume of 

the system and the second one allows one to find the electric charge potential \iq. The 




FIG. 2: The K + /tt + excitation function and energy dependence of the total baryon charge B, 
electrical chemical potential fiQ and volume (sphere radius) for the net strangeness of the system 
5 = 0. Lines are statistical model calculations for classical statistics and s ma x = 1 (i-e. multi- 
strange hyperons are ignored) at two assumptions : (i) dashed lines - fiQ = and correlated sphere 
radius R = 7 fm (as in Fig.l); (ii) solid lines - fj,Q and R are found from the additional conservation 
of baryon and electric charges of the system, Eqs.(@3J and l[4F)[L Experimental points are the same 
as in Fig.l. 
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kaon excitation function recalculated for new obtained volume V = V c is presented in Fig. 2 
together with some thermodynamic quantities. These results are obtained for B =364 and 
Q =142 due to 5% centrality selection of data. 

It is seen that the recalculated K + /tt + ratio does not differ noticeably from that in the 
case of free correlation volume V c , except maybe the SIS energies due to non-vanishing 
fig. The obtained electric chemical potential is not large and it is getting even smaller 
with increasing the invariant colliding energy s 1//2 . However, the baryon conservation is 
fulfilled only for the energy dependent volume. It is noteworthy that the system volume 
behaves non-monotonically having a minimum near s 1 / 2 ~ 10 GeV. This finding may be an 
evidence of the softest point effect at the passage of the phase boundary. In addition, this 
behavior rou ghly correlates with the freeze-out pion volume measured by HBT interference 
experiments [23]. 

Coming back to Fig.l, one can compare the K + /ix + excitation function to those calculated 
within the Maxwell-Boltzmann ensemble but for different values of net strangeness S of the 
hadronic system. It is not of purely academic interest. As was first noted in |24|, in the Gibbs 
phase coexistence of a baryon-rich quark-gluon plasma there is strangeness separation in 
phases, the so-called strangeness distillation effect: While the net strangeness of the system 
is zero, the strangeness of its each component can differ from zero. In particular, in the 
quark subsystem, being a small admixture to the hadronic one, a number of strange quarks 
can be larger than that of strange antiquarks (i.e. S < 0) which should be compensated by 
hadronic phase (i.e. S > 0). A similar situation may occur below the critical temperature 
for systems with crossover deconfinement phase transition. The distillation effect is even 
stronger for lower colliding energies, as recently shown in |2l|. In contrast to the first order 
phase transition, here a small quark admixture can survive during the expansion stage until 
freeze-out. This may give rise to an observable effect if particles from each phase get on 
shell at freeze-out preserving their strangeness. One should remember that recent lattice 
QCD calculations with physical quark masses (3| tell us that at /is = 360 ± 40 MeV and 
T = 162 ± 2 MeV there is the critical end-point where the crossover changes to a phase 
transition of the first order. It means that proceeding from high to lower energies (see 
Fig.l), the hadronic subsystem follows along the curve for K + /n + with S > strangeness 
components, then stops near the end-point and jumps down on the curve with 5 = 0. The 
given parameters of the critical end-point roughly correspond to the bombarding energy 
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~ 30 GeV [2l|| and correlate with the position of the maximum for the relative strangeness 
yield in Fig.l. Therefore, in this scenario the observed maximum is treated as a manifestation 
of the critical end-point rather than a jump due to the latent heat at the first order phase 
transition. A change of the nature of phase transition at the critical end-point may manifest 
itself in other observables as well. Since at this change the dominating hadronic phase of the 
quark-hadron mixture is enriched with strange antiquarks, effects should be more noticeable 
in the K + component. 



V. CONCLUSIONS 



The derivation of the Maxwell-Boltzmann partition function and ensemble averages with 
the exact strangeness conservation has been reconsidered in a general case of non-vanishing 
net strangeness. A new method of recurrence relations was applied to solve the partition 
function. This method has certain advantages as compared to the numerical method used 
for the canonical description of baryon charge conservation of systems with the total 
baryon charge up to about 100. The recurrence relation technique can be successfully applied 
to larger systems, as it was earlier demonstrated by a canonical treatment of statistical 
multifragmentation of excited heavy nuclei 

Using the recurrence relations method we first generalize the calculation scheme for the 
strangeness canonical ensemble to the case of quantum ensemble. Taking into account 
quantum statistics, the K + /n + ratio turns out to be slightly lower at collision energies 
> 6 GeV. 

In the present approach the general behavior of the K + /n + experimental excitation 
function is reproduced roughly to the same level of agreement as that in earlier statistical 
treatments. Uncertainties related to a possible contribution of weak decays and energy 
dependence of dynamic correlation volume are pointed out. While the initial state of the 
equilibrium system is taken from the chemical freeze-out T — /ig curve, an additional account 
of the baryonic and electric charge conservation does not change noticeably the particle 
ratio but provides the thermodynamic consistency of the strangeness canonical approach. 
The obtained unique system volume is energy dependent and exhibits clear minimum near 
s 1 / 2 ~ 10 GeV, which deserves further investigations. 

Nevertheless, the anomalous increase of this ratio near the bombarding energy of about 30 
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GeV, having been interpreted as a sign of deconfinement phase transition, is not reproduced 
by this model. A new scenario based on the strangeness distillation idea has been proposed 
where this anomalous strangeness enhancement is considered as a sign of the evolving system 
passing in vicinity of the critical end-point. Arguments in favor of this hypothesis are given. 
In this respect, multiplicity fluctuations of rare strange particles can be a measure of chemical 
equilibrium [26], and therefore are of special interest. A more convincing test needs delicate 
dynamic calculations which are in progress now. 
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APPENDIX A: RECURSIVE RELATIONS IN CLASSICAL STATISTICS 



Let us prove the recurrence relations (|22|) for the classical perfect gas partition function 
with positive and negative components of the net strangeness. It is easy to see that for the 
given species a in a quantum state p the sum over occupation numbers v av satisfies the 
following identity jicl | 

OO 

E 5 (E s <*>"a> P ' ~ s ) "*P — r e~^ £ap = 

u ap =0 a',p' Va V- 

oo 1 

= e-*°* £ 5(£ s a , u a , p , -(S- s a ))— e"^- . (Al) 

u ap =0 a>,p> Va V 



Taking into account the Kronecker symbol, Eq. pTJ) can be rewritten as 

1 

S: 



4? = 4- E EM E <K E J2KKp-s t )x 



T o'GM T p> {v ap }a&M^ ct6M T P 
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x 1 , e"^^ E » W . (A2) 
Using identity (|A1|) . the summation over i/ ap gives 

4? = }eE Me-^v E «(EE - ( 5 t - M)) x 

x e "^»^ T Ep"^ ( A3) 

n n^i 1 ; 

a£M T p 



By means of Eq. JTJ) the recurrence relations (|A3J) can be obtained directly p, |l 



4? = 4- £ EWeM^i- ( A4 ) 



Thus, the formulae (1221) has been proved. 

Mean occupation numbers are evidently given as 

1 ^ _ ^ _ 1 



Ms = tt E ^(E s ^p - 5 ) "«V Tr \ e-^°* Va * ea * . (A5) 



Z S s v t a.p 11 ^ 

To sum here over iVp' the identity (|Alj) should be used : 



Ms = e-" £ »vi ^ 5(E^^ P " " *a0) e^-^- = (A6) 

Zs <** i], ^ 

= e-^-v ^?si , (A7) 
As 

where Eq.fJTJ has been used for getting the latter equality. 

APPENDIX B: RECURSIVE RELATIONS IN QUANTUM STATISTICS 

Here we will prove the recurrence relations (|3*3~jl for the quantum partition functions with 
positive and negative components of net strangeness. It is readily seen that the following 
identity for occupation number summing is valid 

E u *p e-^ £ap 5(J2 M "a' P > - s±) = 

v ap =0 a',p' 

[S±/\s a \] n>m([S±/\s a \]-l,K a ) 

= E Vaie-^" E e~^^ S{J2\s a/ \u a y-(S ± -l\s a \)) (Bl) 

1=1 u ap =0 a',p' 
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Using the Kronecker symbol, we can rewrite Eq. (|32|) : 

:t) _ i 

? * ~ S: 



4? = }EE M E *( E E I Sq I *V - 5 T ) *Vp< £ 1B2) 



=F a'eA/ T p' {^„} ae M T adM T P 



Applying identity (JBlj) for summing over z/ Q „, we have 



= 4- E EM" E (B3) 

D T a 'eM T p' (=1 

x E *( E E kl *v - (s T - iM)) e~^^ ^>- £ - . 

The use of (EHI) in Eq. (IB3|) immediately gets the recurrence relations licl| 

1 [5t/M1 

4? = ^EEW E ^ te - 4^ w , (b 4 ) 

°=F a£M T P 1=1 

which completes the derivation of Eq. (|33|) . 

The partition function with zero strangeness component has the form 

{fctp}aeM 

This partition function (|B5|) can be represented as the product of sums for every species a 

and then, accounting for the Kronecker symbols, be reduced to the equation 

00 „ 
Z(0) = II E E 6(£ l v e9 -n a )e- f, E> v °*» , (B6) 

oeM n a =0{v ap }> p 

where {v ap }' is the set of occupation numbers for a fixed species a. Formally, sums over 
occupation numbers {z/ ap }' in Eq. (|B6jl are quantum single-particle partition function n a of 
identical particles in the canonical ensemble (see P, [l3|). So, Eq. (|B6|) can be represented as 
follows 

oo n a n a I yl \N a i 

Z {0) = n E E *(£ IN* - na) n J^Tl ■ (B7) 

aeM n a =0 {N al }' 1=1 1=1 1 ah 

After simple transformations we finally get 



[S=/\* a >\] 



oo 

-1 r,l 



Z<®= I] IF a ■ (B8) 

qGM 1=1 

Mean occupation number are given by the following expression 

Ms = 7^ E <KE s a u ap - S) u ap e -^E-*^ = -^^^ • (B9) 
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With (fT9|) . from here we get for mean occupation numbers with a G M T 
Ms — EE 6{-S- + S+-S) z£> x 

s+=0S-=0 

x E 5 ( E E M *V - 5 T ) z> e"^-^ Q G M T (B10) 

Using identity (|BH) and following the derivation of Eqs. (IB2j) - (|B4jl from Eq. (|B10J) . we arrive 
at Eq.JSHl) for the mean occupation numbers. 

In the case of a £ M , the mean occupation numbers are given by the equation 

(v )c = — V v e~ P £ P = -- 9 ln Zi ° aGMa (Bll) 

{ap)s Z(°). P P de ap ' a ^ M °- 

If the partition function (|B8|) for the zero strangeness component is substituted into differ- 
ential relation (|B11|I . then we arrive at Eq. (j39|) for mean occupation numbers with a G Mq. 



